Spatiotemporal genotype replacement of H5N8 avian influenza viruses contributed to H5N1 emergence in 2021/2022 panzootic.

Since 2020, clade 2.3.4.4b highly pathogenic avian influenza H5N8 and H5N1 viruses have swept through continents, posing serious threats to the world. Through comprehensive analyses of epidemiological, genetic, and bird migration data, we found that the dominant genotype replacement of the H5N8 viruses in 2020 contributed to the H5N1 outbreak in the 2021/2022 wave. The 2020 outbreak of the H5N8 G1 genotype instead of the G0 genotype produced reassortment opportunities and led to the emergence of a new H5N1 virus with G1's HA and MP genes. Despite extensive reassortments in the 2021/2022 wave, the H5N1 virus retained the HA and MP genes, causing a significant outbreak in Europe and North America. Furtherly, through the wild bird migration flyways investigation, we found that the temporal-spatial coincidence between the outbreak of the H5N8 G1 virus and the bird autumn migration may have expanded the H5 viral spread, which may be one of the main drivers of the emergence of the 2020-2022 H5 panzootic.IMPORTANCESince 2020, highly pathogenic avian influenza (HPAI) H5 subtype variants of clade 2.3.4.4b have spread across continents, posing unprecedented threats globally. However, the factors promoting the genesis and spread of H5 HPAI viruses remain unclear. Here, we found that the spatiotemporal genotype replacement of H5N8 HPAI viruses contributed to the emergence of the H5N1 variant that caused the 2021/2022 panzootic, and the viral evolution in poultry of Egypt and surrounding area and autumn bird migration from the Russia-Kazakhstan region to Europe are important drivers of the emergence of the 2020-2022 H5 panzootic. These findings provide important targets for early warning and could help control the current and future HPAI epidemics.

DiffDomain enables identification of structurally reorganized topologically associating domains.

Topologically associating domains (TADs) are critical structural units in three-dimensional genome organization of mammalian genome. Dynamic reorganizations of TADs between health and disease states are associated with essential genome functions. However, computational methods for identifying reorganized TADs are still in the early stages of development. Here, we present DiffDomain, an algorithm leveraging high-dimensional random matrix theory to identify structurally reorganized TADs using high-throughput chromosome conformation capture (Hi-C) contact maps. Method comparison using multiple real Hi-C datasets reveals that DiffDomain outperforms alternative methods for false positive rates, true positive rates, and identifying a new subtype of reorganized TADs. Applying DiffDomain to Hi-C data from different cell types and disease states demonstrates its biological relevance. Identified reorganized TADs are associated with structural variations and epigenomic changes such as changes in CTCF binding sites. By applying to a single-cell Hi-C data from mouse neuronal development, DiffDomain can identify reorganized TADs between cell types with reasonable reproducibility using pseudo-bulk Hi-C data from as few as 100 cells per condition. Moreover, DiffDomain reveals differential cell-to-population variability and heterogeneous cell-to-cell variability in TADs. Therefore, DiffDomain is a statistically sound method for better comparative analysis of TADs using both Hi-C and single-cell Hi-C data.

Genetic determinants of IgG antibody response to COVID-19 vaccination.

Human humoral immune responses to SARS-CoV-2 vaccines exhibit substantial inter-individual variability and have been linked to vaccine efficacy. To elucidate the underlying mechanism behind this variability, we conducted a genome-wide association study (GWAS) on the anti-spike IgG serostatus of UK Biobank participants who were previously uninfected by SARS-CoV-2 and had received either the first dose (n = 54,066) or the second dose (n = 46,232) of COVID-19 vaccines. Our analysis revealed significant genome-wide associations between the IgG antibody serostatus following the initial vaccine and human leukocyte antigen (HLA) class II alleles. Specifically, the HLA-DRB1∗13:02 allele (MAF = 4.0%, OR = 0.75, p = 2.34e-16) demonstrated the most statistically significant protective effect against IgG seronegativity. This protective effect was driven by an alteration from arginine (Arg) to glutamic acid (Glu) at position 71 on HLA-DRß1 (p = 1.88e-25), leading to a change in the electrostatic potential of pocket 4 of the peptide binding groove. Notably, the impact of HLA alleles on IgG responses was cell type specific, and we observed a shared genetic predisposition between IgG status and susceptibility/severity of COVID-19. These results were replicated within independent cohorts where IgG serostatus was assayed by two different antibody serology tests. Our findings provide insights into the biological mechanism underlying individual variation in responses to COVID-19 vaccines and highlight the need to consider the influence of constitutive genetics when designing vaccination strategies for optimizing protection and control of infectious disease across diverse populations.

Investigating linear and nonlinear associations of LDL cholesterol with incident chronic kidney disease, atherosclerotic cardiovascular disease and all-cause mortality: A prospective and Mendelian randomization study.

BACKGROUND AND AIMS: Observational studies suggest potential nonlinear associations of low-density lipoprotein cholesterol (LDL-C) with cardio-renal diseases and mortality, but the causal nature of these associations is unclear. We aimed to determine the shape of causal relationships of LDL-C with incident chronic kidney disease (CKD), atherosclerotic cardiovascular disease (ASCVD) and all-cause mortality, and to evaluate the absolute risk of adverse outcomes contributed by LDL-C itself. METHODS: Observational analysis and one-sample Mendelian randomization (MR) with linear and nonlinear assumptions were performed using the UK Biobank of >0.3 million participants with no reported prescription of lipid-lowering drugs. Two-sample MR on summary-level data from the Global Lipid Genetics Consortium (N = 296,680) and the CKDGen (N = 625,219) was employed to replicate the relationship for kidney traits. The 10-year probabilities of the outcomes was estimated by integrating the MR and Cox models. RESULTS: Observationally, participants with low LDL-C were significantly associated with a decreased risk of ASCVD, but an increased risk of CKD and all-cause mortality. Univariable MR showed an inverse total effect of LDL-C on incident CKD (HR [95% CI]:0.84 [0.73-0.96]; p = 0.011), a positive effect on ASCVD (1.41 [1.29-1.53]; p<0.001), and no significant causal effect on all-cause mortality. Multivariable MR, controlling for high-density lipoprotein cholesterol (HDL-C) and triglycerides, identified a positive direct effect on ASCVD (1.32 [1.18-1.47]; p<0.001), but not on CKD and all-cause mortality. These results indicated that genetically predicted low LDL-C had an inverse indirect effect on CKD mediated by HDL-C and triglycerides, which was validated by a two-sample MR analysis using summary-level data from the Global Lipid Genetics Consortium (N = 296,680) and the CKDGen consortium (N = 625,219). Suggestive evidence of a nonlinear causal association between LDL-C and CKD was found. The 10-year probability curve showed that LDL-C concentrations below 3.5 mmol/L were associated with an increased risk of CKD. CONCLUSIONS: In the general population, lower LDL-C was causally associated with lower risk of ASCVD, but appeared to have a trade-off for an increased risk of CKD, with not much effect on all-cause mortality. LDL-C concentration below 3.5 mmol/L may increase the risk of CKD.

Effect of seasonal coronavirus immune imprinting on the immunogenicity of inactivated COVID-19 vaccination.

Background: Pre-existing cross-reactive immunity among different coronaviruses, also termed immune imprinting, may have a comprehensive impact on subsequent SARS-CoV-2 infection and COVID-19 vaccination effectiveness. Here, we aim to explore the interplay between pre-existing seasonal coronaviruses (sCoVs) antibodies and the humoral immunity induced by COVID-19 vaccination. Methods: We first collected serum samples from healthy donors prior to COVID-19 pandemic and individuals who had received COVID-19 vaccination post-pandemic in China, and the levels of IgG antibodies against sCoVs and SARS-CoV-2 were detected by ELISA. Wilcoxon rank sum test and chi-square test were used to compare the difference in magnitude and seropositivity rate between two groups. Then, we recruited a longitudinal cohort to collect serum samples before and after COVID-19 vaccination. The levels of IgG antibodies against SARS-CoV-2 S, S1, S2 and N antigen were monitored. Association between pre-existing sCoVs antibody and COVID-19 vaccination-induced antibodies were analyzed by Spearman rank correlation. Results: 96.0% samples (339/353) showed the presence of IgG antibodies against at least one subtype of sCoVs. 229E and OC43 exhibited the highest seroprevalence rates at 78.5% and 72.0%, respectively, followed by NL63 (60.9%) and HKU1 (52.4%). The levels of IgG antibodies against two ß coronaviruses (OC43 and HKU1) were significantly higher in these donors who had inoculated with COVID-19 vaccines compared to pre-pandemic healthy donors. However, we found that COVID-19 vaccine-induced antibody levels were not significant different between two groups with high levelor low level of pre-existing sCoVs antibody among the longitudinal cohort. Conclusion: We found a high prevalence of antibodies against sCoVs in Chinese population. The immune imprinting by sCoVs could be reactivated by COVID-19 vaccination, but it did not appear to be a major factor affecting the immunogenicity of COVID-19 vaccine. These findings will provide insights into understanding the impact of immune imprinting on subsequent multiple shots of COVID-19 vaccines.

Characterizing the dynamics of BCR repertoire from repeated influenza vaccination.

Yearly epidemics of seasonal influenza cause an enormous disease burden around the globe. An understanding of the rules behind the immune response with repeated vaccination still presents a significant challenge, which would be helpful for optimizing the vaccination strategy. In this study, 34 healthy volunteers with 16 vaccinated were recruited, and the dynamics of the BCR repertoire for consecutive vaccinations in two seasons were tracked. In terms of diversity, length, network, V and J gene segments usage, somatic hypermutation (SHM) rate and isotype, it was found that the overall changes were stronger in the acute phase of the first vaccination than the second vaccination. However, the V gene segments of IGHV4-39, IGHV3-9, IGHV3-7 and IGHV1-69 were amplified in the acute phase of the first vaccination, with IGHV3-7 dominant. On the other hand, for the second vaccination, the changes were dominated by IGHV1-69, with potential for coding broad neutralizing antibody. Additional analysis indicates that the application of V gene segment for IGHV3-7 in the acute phase of the first vaccination was due to the elevated usage of isotypes IgM and IgG3. While for IGHV1-69 in the second vaccination, it was contributed by isotypes IgG1 and IgG2. Finally, 41 public BCR clusters were identified in the vaccine group, with both IGHV3-7 and IGHV1-69 were involved and representative complementarity determining region 3 (CDR3) motifs were characterized. This study provides insights into the immune response dynamics following repeated influenza vaccination in humans and can inform universal vaccine design and vaccine strategies in the future.

Effect of Rapid Urbanization in Mainland China on the Seasonal Influenza Epidemic: Spatiotemporal Analysis of Surveillance Data From 2010 to 2017.

BACKGROUND: The world is undergoing an unprecedented wave of urbanization. However, the effect of rapid urbanization during the early or middle stages of urbanization on seasonal influenza transmission remains unknown. Since about 70% of the world population live in low-income countries, exploring the impact of urbanization on influenza transmission in urbanized countries is significant for global infection prediction and prevention. OBJECTIVE: The aim of this study was to explore the effect of rapid urbanization on influenza transmission in China. METHODS: We performed spatiotemporal analyses of province-level influenza surveillance data collected in Mainland China from April 1, 2010, to March 31, 2017. An agent-based model based on hourly human contact-related behaviors was built to simulate the influenza transmission dynamics and to explore the potential mechanism of the impact of urbanization on influenza transmission. RESULTS: We observed persistent differences in the influenza epidemic attack rates among the provinces of Mainland China across the 7-year study period, and the attack rate in the winter waves exhibited a U-shaped relationship with the urbanization rates, with a turning point at 50%-60% urbanization across Mainland China. Rapid Chinese urbanization has led to increases in the urban population density and percentage of the workforce but decreases in household size and the percentage of student population. The net effect of increased influenza transmission in the community and workplaces but decreased transmission in households and schools yielded the observed U-shaped relationship. CONCLUSIONS: Our results highlight the complicated effects of urbanization on the seasonal influenza epidemic in China. As the current urbanization rate in China is approximately 59%, further urbanization with no relevant interventions suggests a worrisome increasing future trend in the influenza epidemic attack rate.

Prediction of red blood cell transfusion after orthopedic surgery using an interpretable machine learning framework.

Objective: Postoperative red blood cell (RBC) transfusion is widely used during the perioperative period but is often associated with a high risk of infection and complications. However, prediction models for RBC transfusion in patients with orthopedic surgery have not yet been developed. We aimed to identify predictors and constructed prediction models for RBC transfusion after orthopedic surgery using interpretable machine learning algorithms. Methods: This retrospective cohort study reviewed a total of 59,605 patients undergoing orthopedic surgery from June 2013 to January 2019 across 7 tertiary hospitals in China. Patients were randomly split into training (80%) and test subsets (20%). The feature selection method of recursive feature elimination (RFE) was used to identify an optimal feature subset from thirty preoperative variables, and six machine learning algorithms were applied to develop prediction models. The Shapley Additive exPlanations (SHAP) value was employed to evaluate the contribution of each predictor towards the prediction of postoperative RBC transfusion. For simplicity of the clinical utility, a risk score system was further established using the top risk factors identified by machine learning models. Results: Of the 59,605 patients with orthopedic surgery, 19,921 (33.40%) underwent postoperative RBC transfusion. The CatBoost model exhibited an AUC of 0.831 (95% CI: 0.824-0.836) on the test subset, which significantly outperformed five other prediction models. The risk of RBC transfusion was associated with old age (>60 years) and low RBC count (<4.0 × 1012/L) with clear threshold effects. Extremes of BMI, low albumin, prolonged activated partial thromboplastin time, repair and plastic operations on joint structures were additional top predictors for RBC transfusion. The risk score system derived from six risk factors performed well with an AUC of 0.801 (95% CI: 0.794-0.807) on the test subset. Conclusion: By applying an interpretable machine learning framework in a large-scale multicenter retrospective cohort, we identified novel modifiable risk factors and developed prediction models with good performance for postoperative RBC transfusion in patients undergoing orthopedic surgery. Our findings may allow more precise identification of high-risk patients for optimal control of risk factors and achieve personalized RBC transfusion for orthopedic patients.

Variation in synonymous evolutionary rates in the SARS-CoV-2 genome.

Introduction: Coronavirus disease 2019 is an infectious disease caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Influential variants and mutants of this virus continue to emerge, and more effective virus-related information is urgently required for identifying and predicting new mutants. According to earlier reports, synonymous substitutions were considered phenotypically silent; thus, such mutations were frequently ignored in studies of viral mutations because they did not directly cause amino acid changes. However, recent studies have shown that synonymous substitutions are not completely silent, and their patterns and potential functional correlations should thus be delineated for better control of the pandemic. Methods: In this study, we estimated the synonymous evolutionary rate (SER) across the SARS-CoV-2 genome and used it to infer the relationship between the viral RNA and host protein. We also assessed the patterns of characteristic mutations found in different viral lineages. Results: We found that the SER varies across the genome and that the variation is primarily influenced by codon-related factors. Moreover, the conserved motifs identified based on the SER were found to be related to host RNA transport and regulation. Importantly, the majority of the existing fixed-characteristic mutations for five important virus lineages (Alpha, Beta, Gamma, Delta, and Omicron) were significantly enriched in partially constrained regions. Discussion: Taken together, our results provide unique information on the evolutionary and functional dynamics of SARS-CoV-2 based on synonymous mutations and offer potentially useful information for better control of the SARS-CoV-2 pandemic.

Mechanisms for the circulation of influenza A(H3N2) in China: A spatiotemporal modelling study.

Circulation of seasonal influenza is the product of complex interplay among multiple drivers, yet characterizing the underlying mechanism remains challenging. Leveraging the diverse seasonality of A(H3N2) virus and abundant climatic space across regions in China, we quantitatively investigated the relative importance of population susceptibility, climatic factors, and antigenic change on the dynamics of influenza A(H3N2) through an integrative modelling framework. Specifically, an absolute humidity driven multiscale transmission model was constructed for the 2013/2014, 2014/2015 and 2016/2017 influenza seasons that were dominated by influenza A(H3N2). We revealed the variable impact of absolute humidity on influenza transmission and differences in the occurring timing and magnitude of antigenic change for those three seasons. Overall, the initial population susceptibility, climatic factors, and antigenic change explained nearly 55% of variations in the dynamics of influenza A(H3N2). Specifically, the additional variation explained by the initial population susceptibility, climatic factors, and antigenic change were at 33%, 26%, and 48%, respectively. The vaccination program alone failed to fully eliminate the summer epidemics of influenza A(H3N2) and non-pharmacological interventions were needed to suppress the summer circulation. The quantitative understanding of the interplay among driving factors on the circulation of influenza A(H3N2) highlights the importance of simultaneous monitoring of fluctuations for related factors, which is crucial for precise and targeted prevention and control of seasonal influenza.

A network view of human immune system and virus-human interaction.

The immune system is highly networked and complex, which is continuously changing as encountering old and new pathogens. However, reductionism-based researches do not give a systematic understanding of the molecular mechanism of the immune response and viral pathogenesis. Here, we present HUMPPI-2022, a high-quality human protein-protein interaction (PPI) network, containing > 11,000 protein-coding genes with > 78,000 interactions. The network topology and functional characteristics analyses of the immune-related genes (IRGs) reveal that IRGs are mostly located in the center of the network and link genes of diverse biological processes, which may reflect the gene pleiotropy phenomenon. Moreover, the virus-human interactions reveal that pan-viral targets are mostly hubs, located in the center of the network and enriched in fundamental biological processes, but not for coronavirus. Finally, gene age effect was analyzed from the view of the host network for IRGs and virally-targeted genes (VTGs) during evolution, with IRGs gradually became hubs and integrated into host network through bridging functionally differentiated modules. Briefly, HUMPPI-2022 serves as a valuable resource for gaining a better understanding of the composition and evolution of human immune system, as well as the pathogenesis of viruses.

Complex q-rung orthopair fuzzy Frank aggregation operators and their application to multi-attribute decision making.

The complex q-rung orthopair fuzzy sets (Cq-ROFSs) can serve as a generalization of q-rung orthopair fuzzy sets (q-ROFSs) and complex fuzzy sets FS (CFSs). Cq-ROFSs provide more freedom for people handling uncertainty and vagueness by the truth and falsity grades on the condition that the sum of the q-powers of the real part and imaginary part is within the unit interval. Further, Frank operational laws are an extended form of Archimedes' T mode and Archimedes' S mode and Frank aggregation operators have a certain parameter which makes them more flexible and more generalized than many other aggregation operators in the process of information fusion. The objectives of this paper are to extend the Frank operations to the complex q-rung orthopair fuzzy environment and to introduce their score function and accuracy function. Meanwhile, some complex q-rung fuzzy Frank aggregation operators are developed, such as the complex q-rung orthopair fuzzy Frank weighted averaging (Cq-ROFFWA) operator, the complex q-rung orthopair fuzzy Frank weighted geometric (Cq-ROFFWG) operator, and the complex q-rung orthopair fuzzy Frank ordered weighted averaging (Cq-ROFFOWA) operator, and their special cases are discussed. In addition, an innovative MADM method is introduced according to the propounded operators to deal with multi-attribute decision-making problems under the complex q-rung orthopair fuzzy environment. Consequently, the practicability and effectiveness of the created methods are proposed by parameter exploration and comparative analysis.

Transmission Patterns of Seasonal Influenza in China between 2010 and 2018.

Background Understanding the transmission source, pattern, and mechanism of infectious diseases is essential for targeted prevention and control. Though it has been studied for many years, the detailed transmission patterns and drivers for the seasonal influenza epidemics in China remain elusive. Methods In this study, utilizing a suite of epidemiological and genetic approaches, we analyzed the updated province-level weekly influenza surveillance, sequence, climate, and demographic data between 1 April 2010 and 31 March 2018 from continental China, to characterize detailed transmission patterns and explore the potential initiating region and drivers of the seasonal influenza epidemics in China. Results An annual cycle for influenza A(H1N1)pdm09 and B and a semi-annual cycle for influenza A(H3N2) were confirmed. Overall, the seasonal influenza A(H3N2) virus caused more infection in China and dominated the summer season in the south. The summer season epidemics in southern China were likely initiated in the "Lingnan" region, which includes the three most southern provinces of Hainan, Guangxi, and Guangdong. Additionally, the regions in the south play more important seeding roles in maintaining the circulation of seasonal influenza in China. Though intense human mobility plays a role in the province-level transmission of influenza epidemics on a temporal scale, climate factors drive the spread of influenza epidemics on both the spatial and temporal scales. Conclusion The surveillance of seasonal influenza in the south, especially the "Lingnan" region in the summer, should be strengthened. More broadly, both the socioeconomic and climate factors contribute to the transmission of seasonal influenza in China. The patterns and mechanisms revealed in this study shed light on the precise forecasting, prevention, and control of seasonal influenza in China and worldwide.

Subtypes specified environmental dependence of seasonal influenza virus.

Understanding the role of meteorological factors in the transmission dynamics of respiratory infectious diseases remains challenging. Our study was to comprehensively investigate the nonlinear effects of environmental factors on influenza transmission, based on multi-region surveillance data from mainland China. An approach related to time-varying reproduction number (Rt) was proposed, which extracts the environment-related components from Rt to estimate the relationship between environmental factors and influenza transmission based on a mixed-effects regression model. Nonlinear relationships for absolute humidity (the lowest transmission was observed at absolute humidity of 12 g/m3) and mean temperature (the lowest transmission was observed at the mean temperature of 18 °C) with influenza transmission were observed. Influenza transmission holds almost constant with the average precipitation below 1 mm or sunshine hour below 9 h/day, but increases for the precipitation and decreases for the sunshine hour afterward. The environmental dependence varies across subtypes: A(H3N2) maintains relatively higher transmission in high temperature and humidity conditions, compared with other influenza subtypes. Overall, the subtypes specified environmental dependence of influenza transmission could explain 23.1 %, 29.2 % and 27.1 % of the variations for A(H1N1)pdm09, A(H3N2) and B-lineage in China. The projected seasonal transmission rates based on our approach could be used as a valuable seasonal proxy to model the influenza dynamics under various meteorological spaces. Finally, our approach is also applicable to obtain novel insights into the impact of environmental factors on other respiratory infectious diseases.

Vaccination strategy for preventing the spread of SARS-CoV-2 in the limited supply condition: A mathematical modeling study.

To mitigate SARS-CoV-2 transmission, vaccines have been urgently approved. With their limited availability, it is critical to distribute the vaccines reasonably. We simulated the SARS-CoV-2 transmission for 365 days over four intervention periods: free transmission, structural mitigation, personal mitigation, and vaccination. Sensitivity analyses were performed to obtain robust results. We further evaluated two proposed vaccination allocations, including one-dose-high-coverage and two-doses-low-coverage, when the supply was low. 33.35% (infection rate, 2.68 in 10 million people) and 40.54% (2.36) of confirmed cases could be avoided as the nonpharmaceutical interventions (NPIs) adherence rate rose from 50% to 70%. As the vaccination coverage reached 60% and 80%, the total infections could be reduced by 32.72% and 41.19%, compared to the number without vaccination. When the durations of immunity were 90 and 120 days, the infection rates were 2.67 and 2.38. As the asymptomatic infection rate rose from 30% to 50%, the infection rate increased 0.92 (SD, 0.16) times. Conditioned on 70% adherence rate, with the same amount of limited available vaccines, the 20% and 40% vaccination coverage of one-dose-high-coverage, the infection rates were 2.70 and 2.35; corresponding to the two-doses-low-coverage with 10% and 20% vaccination coverage, the infection rates were 3.22 and 2.92. Our results indicated as the duration of immunity prolonged, the second wave of SARS-CoV-2 would be delayed and the scale would be declined. On average, the total infections in two-doses-low-coverage was 1.48 times (SD, 0.24) as high as that in one-dose-high-coverage. It is crucial to encourage people in order to improve vaccination coverage and establish immune barriers. Particularly when the supply is limited, a wiser strategy to prevent SARS-CoV-2 is equally distributing doses to the same number of individuals. Besides vaccination, NPIs are equally critical to the prevention of widespread of SARS-CoV-2.

The effectiveness of governmental nonpharmaceutical interventions against COVID-19 at controlling seasonal influenza transmission: an ecological study.

BACKGROUND: A range of strict nonpharmaceutical interventions (NPIs) were implemented in many countries to combat the coronavirus 2019 (COVID-19) pandemic. These NPIs may also be effective at controlling seasonal influenza virus infections, as influenza viruses have the same transmission path as severe acute respiratory syndrome coronavirus 2. The aim of this study was to evaluate the effects of different NPIs on the control of seasonal influenza. METHODS: Data for 14 NPIs implemented in 33 countries and the corresponding influenza virological surveillance data were collected. The influenza suppression index was calculated as the difference between the influenza positivity rate during its period of decline from 2019 to 2020 and during the influenza epidemic seasons in the previous 9 years. A machine learning model was developed using an extreme gradient boosting tree regressor to fit the NPI and influenza suppression index data. The SHapley Additive exPlanations tool was used to characterize the NPIs that suppressed the transmission of influenza. RESULTS: Of all NPIs tested, gathering limitations had the greatest contribution (37.60%) to suppressing influenza transmission during the 2019-2020 influenza season. The three most effective NPIs were gathering limitations, international travel restrictions, and school closures. For these three NPIs, their intensity threshold required to generate an effect were restrictions on the size of gatherings less than 1000 people, ban of travel to all regions or total border closures, and closing only some categories of schools, respectively. There was a strong positive interaction effect between mask-wearing requirements and gathering limitations, whereas merely implementing a mask-wearing requirement, and not other NPIs, diluted the effectiveness of mask-wearing requirements at suppressing influenza transmission. CONCLUSIONS: Gathering limitations, ban of travel to all regions or total border closures, and closing some levels of schools were found to be the most effective NPIs at suppressing influenza transmission. It is recommended that the mask-wearing requirement be combined with gathering limitations and other NPIs. Our findings could facilitate the precise control of future influenza epidemics and other potential pandemics.

Drivers and forecasts of multiple waves of the coronavirus disease 2019 pandemic: A systematic analysis based on an interpretable machine learning framework.

Coronavirus disease 2019 (COVID-19) has become a global pandemic and continues to prevail with multiple rebound waves in many countries. The driving factors for the spread of COVID-19 and their quantitative contributions, especially to rebound waves, are not well studied. Multidimensional time-series data, including policy, travel, medical, socioeconomic, environmental, mutant and vaccine-related data, were collected from 39 countries up to 30 June 2021, and an interpretable machine learning framework (XGBoost model with Shapley Additive explanation interpretation) was used to systematically analyze the effect of multiple factors on the spread of COVID-19, using the daily effective reproduction number as an indicator. Based on a model of the pre-vaccine era, policy-related factors were shown to be the main drivers of the spread of COVID-19, with a contribution of 60.81%. In the post-vaccine era, the contribution of policy-related factors decreased to 28.34%, accompanied by an increase in the contribution of travel-related factors, such as domestic flights, and contributions emerged for mutant-related (16.49%) and vaccine-related (7.06%) factors. For single-peak countries, the dominant ones were policy-related factors during both the rising and fading stages, with overall contributions of 33.7% and 37.7%, respectively. For double-peak countries, factors from the rebound stage contributed 45.8% and policy-related factors showed the greatest contribution in both the rebound (32.6%) and fading (25.0%) stages. For multiple-peak countries, the Delta variant, domestic flights (current month) and the daily vaccination population are the three greatest contributors (8.12%, 7.59% and 7.26%, respectively). Forecasting models to predict the rebound risk were built based on these findings, with accuracies of 0.78 and 0.81 for the pre- and post-vaccine eras, respectively. These findings quantitatively demonstrate the systematic drivers of the spread of COVID-19, and the framework proposed in this study will facilitate the targeted prevention and control of the ongoing COVID-19 pandemic.

Susceptibility identification for seasonal influenza A/H3N2 based on baseline blood transcriptome.

Introduction: Influenza susceptibility difference is a widely existing trait that has great practical significance for the accurate prevention and control of influenza. Methods: Here, we focused on the human susceptibility to the seasonal influenza A/H3N2 of healthy adults at baseline level. Whole blood expression data for influenza A/H3N2 susceptibility from GEO were collected firstly (30 symptomatic and 19 asymptomatic). Then to explore the differences at baseline, a suite of systems biology approaches - the differential expression analysis, co-expression network analysis, and immune cell frequencies analysis were utilized. Results: We found the baseline condition, especially immune condition between symptomatic and asymptomatic, was different. Co-expression module that is positively related to asymptomatic is also related to immune cell type of naïve B cell. Function enrichment analysis showed significantly correlation with "B cell receptor signaling pathway", "immune response-activating cell surface receptor signaling pathway" and so on. Also, modules that are positively related to symptomatic are also correlated to immune cell type of neutrophils, with function enrichment analysis showing significantly correlations with "response to bacterium", "inflammatory response", "cAMP-dependent protein kinase complex" and so on. Responses of symptomatic and asymptomatic hosts after virus exposure show differences on resisting the virus, with more effective frontline defense for asymptomatic hosts. A prediction model was also built based on only baseline transcription information to differentiate symptomatic and asymptomatic population with accuracy of 0.79. Discussion: The results not only improve our understanding of the immune system and influenza susceptibility, but also provide a new direction for precise and targeted prevention and therapy of influenza.

Heavy chain sequence-based classifier for the specificity of human antibodies.

Antibodies specifically bind to antigens and are an essential part of the immune system. Hence, antibodies are powerful tools in research and diagnostics. High-throughput sequencing technologies have promoted comprehensive profiling of the immune repertoire, which has resulted in large amounts of antibody sequences that remain to be further analyzed. In this study, antibodies were downloaded from IMGT/LIGM-DB and Sequence Read Archive databases. Contributing features from antibody heavy chains were formulated as numerical inputs and fed into an ensemble machine learning classifier to classify the antigen specificity of six classes of antibodies, namely anti-HIV-1, anti-influenza virus, anti-pneumococcal polysaccharide, anti-citrullinated protein, anti-tetanus toxoid and anti-hepatitis B virus. The classifier was validated using cross-validation and a testing dataset. The ensemble classifier achieved a macro-average area under the receiver operating characteristic curve (AUC) of 0.9246 from the 10-fold cross-validation, and 0.9264 for the testing dataset. Among the contributing features, the contribution of the complementarity-determining regions was 53.1% and that of framework regions was 46.9%, and the amino acid mutation rates occupied the first and second ranks among the top five contributing features. The classifier and insights provided in this study could promote the mechanistic study, isolation and utilization of potential therapeutic antibodies.

Assessment of the global circulation and endemicity of dengue.

Dengue is a significant public health issue, affecting hundreds of millions of people worldwide. As it is spreading from tropical and subtropical zones, some regions previously recognised as non-endemic are at risk of becoming endemic. However, the global circulation of dengue is not fully understood and quantitative measurements of endemicity levels are lacking, posing an obstacle in the precise control of dengue spread. In this study, a sequence-based pipeline was designed based on random sampling to study the transmission of dengue. The limited intercontinental transmission was identified, while regional circulation of dengue was quantified in terms of importation, local circulation and exportation. Additionally, hypo- and hyper-endemic regions were identified using a new metric, with the former characterised by low local circulation and increased importation, whereas the latter by high local circulation and reduced importation. In this study, the global circulation pattern of dengue was examined and a sequence-based endemicity measurement was proposed, which will be helpful for future surveillance and targeted control of dengue.